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Abstract 

In this paper we investigate the integrabihty of two-dimensional partial difference equations 
using the newly developed techniques of study of the degree of the iterates. We show that 
while for generic, nonintegrable equations, the degree grows exponentially fast, for integrable 
lattice equations the degree growth is polynomial. The growth criterion is used in order to 
obtain the integrable deautonomisations of the equations examined. In the case of linearisable 
lattice equations we show that the degree growth is slower than in the case of equations 
integrable through Inverse Scattering Transform techniques. 



The study of integrability of nonlinear evolution equations has spurred the development of 
efficient tools for its detection. The ARS [1] conjecture was formulated originally for partial 
differential equations and related integrability to the Painleve property. In the discrete 
domain the singularity confinement [2] property was discovered while studying the lattice 
KdV equation and the singularities that can appear spontaneously during the evolution. The 
singularity confinement has been a most useful discrete integrability criterion in the sense 
that it is a necessary condition for lattice equations to be integrable by Inverse Scattering 
Transform (1ST) methods. However, it has turned out [3] that singularity confinement is not 
sufficient for integrability and thus its use as an integrability detector must be subject to 
particular caution. 

Another property of integrable discrete systems, namely the growth of the degree of the 
iterates [4], has, in the long run, proven to be a reliable integrability detector. The main 
idea goes back to Arnold [5] and Veselov [6]. As Veselov summarized it: "integrability has 
an essential correlation with the weak growth of certain characteristics" . The characteristic 
quantity which can be easily obtained and computed for a rational mapping is the degree 
of the numerators or denominators of (the irreducible forms of) the iterates of some initial 
condition. (In order to obtain the degree one must introduce homogeneous coordinates and 
compute the homogeneity degree). Those ideas were refined by Viallet and collaborators 
[7,8], leading to the introduction of the notion of algebraic entropy. The latter is defined as 
E = lim^^oo log((i„)/n where dn is the degree of the n-th iterate. A generic, nonintegrable, 
mapping leads to exponential growth of the degrees of the iterates and thus has a nonzero 
algebraic entropy, while an integrable mapping has zero algebraic entropy. As we have shown 
in a previous work [9] , this is too crude an estimate. The degree growth contains information 
that can be an indication as to the precise integration method to be used and thus should 
be studied in detail. (At this point, we must stress that, as was already pointed out in [8], 
the degrees of the iterates are not invariant under transformation of the variables. However 
the degree growth is invariant and characterises the system at hand). 

In previous works of ours we have applied the techniques of degree growth to the study of 
one-dimensional mappings [9,10]. A first important conclusion of these studies was the con- 
firmation of the singularity confinement results [11] on the derivation of discrete Painleve 
equations. We have shown that, when singularity confinement is used for the deautonomi- 
sation of an integrable autonomous mapping, the condition obtained is identical to the one 
found by requiring nonexponential growth of the degrees of the iterates. (The terms "degrees 
of the iterates" in the above sentence and in the rest of the paper must be understood as 
the common homogeneity degree of the numerators and denominators of their irreducible 
forms, obtained through the introduction of the homogeneous coordinates). This not only 
confirms the results previously obtained for discrete Painleve equations, but also suggests a 
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dual strategy for the study of discrete integrability based on the combined use of singularity 
confinement and study of degree growth. The second result [9] was that mappings which are 
linearisable are associated to a degree growth slower than the ones integrable through 1ST 
techniques. Thus, the detailed study of the degree is not only an indication of integrability 
but also of the integration method. 

In this paper, we apply the techniques of degree growth to two-dimensional partial difference 
equations. We shall show that the main conclusions from the study of one-dimensional 
mappings carry over to the two-dimensional case in a rather straightforward way. 

Let us start with the examination of the equation that serves as a paradigm in all integrability 
studies, namely KdV, the discrete form of which is [12,13]: 
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(Incidentally, this is precisely the equation we have studied in [2], while investigating the 
singularity confinement property.) The study of the degree growth of the iterates in the case 
of a 2-dimensional lattice is substantially more difficult than that of the 1-dimensional case. 
It is thus very important to make the right choices from the outset. Here are the initial 
conditions we choose: on the line m = we take of the form — Pn/q while on the line 
n = we choose Xq" = Tm/q (with ro = Po). We assign to q and the p's, r's the same degree 
of homogeneity. Then we compute the iterates of X using (1) and calculate the degree of 
homogeneity in p, q, r at the various points of the lattice. Here is what we find: 



19 31 41 51 



13 19 25 31 
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At this point we must indicate how the analytical expression for the degree can be obtained. 
First we compute several points on the lattice which allow us to have a good guess at how 
the degree behaves. In the particular case of a 2-dimensional discrete equation relating four 
points on an elementary square like (1), and with the present choice of initial conditions 
(and given our experience on 1-dimensional mappings) we can reasonably surmise that the 
dominant behaviour of the degree will be of the form oc mn. Moreover the subdominant 
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terms must be symmetric in m, n and at most linear. With those indications it is possible to 
"guess" the expression = Amn — 2 max(m, n) + l (for mn ^ 0) and subsequently calculate 
some more points in order to check its validity. This procedure will be used throughout this 
paper. 

So the lattice KdV equation leads, quite expectedly, to a polynomial growth in the degrees 
of the iterates. Let us now turn to the more interesting question of deautonomisation. The 
form (1) of KdV is not very convenient and thus we shall study its potential form [14]: 

m+1 _ m _| fn ln\ 

- -^n ^ m+1 _ „m " ^"^^ 

(The name 'potential' is given here in analogy to the continuous case: the dependent variable 
X of equation (2) is related to the dependent variable X of equation (1) through a;^"*"^— a^^^.^ = 
and (1) is recovered exactly if z^^—l). The deautonomisation we are referring to consists 
in finding an explicit m, n dependence of which is compatible with integrability. Let us 
first compute the degrees of the iterates for constant z: 
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The degree is given simply by = mn + 1. Assuming a generic (m, n) dependence for 
z we obtain the following successive degrees: 
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We remark readily that the degrees form a Pascal triangle i.e. they are identical to the 
binomial coefficients, leading to an exponential growth at least on a strip along the diagonal. 
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The way to obtain an integrable deautonomisation is to require that the degrees obtained in 
the autonomous and nonautonomous cases be identical. The first constraint can be obtained 
by reducing the degree of from 6 to 5. As a matter of fact, starting from the initial 
conditions a;° = p^/q , = r^/g (with tq = po) we obtain x\ = ipiPo-pori-z^q^)/{q{pi - 
ri)), Xi = Qz/{qQ2) where Qk is a polynomial of degree fc, and a similar expression for x\. 
Computing x^ we find X2 — Qe/iqipi — i'i)Q4)- It is impossible for q to divide Qq for generic 
initial conditions. However, requiring {pi — ri) to be a factor of Qq we find the constraint 
z\ — Zq — Zi + Zq = 0. The relation of this result to singularity confinement is quite easy to 
perceive. The singularity corresponding to g = is indeed a fixed singularity: it exists for 
all (n, m)'s where either n or m are equal to zero. On the other hand the singularity related 
to pi — ri =0 appears only at a certain iteration and is thus movable. The fact that with 
the proper choice of the denominator factors out, is precisely what one expects for the 
singularity to be confined. 



m+l _ m+1 _ m , m _ rv (o\ 



Requiring that z satisfy 

suffices to reduce the degrees of all higher x's to those of the autonomous case. The solution 
of (3) is z!^ = f{n) + g(m) where /, g are two arbitrary functions. This form of z^^ is 
precisely the one obtained in the analysis of convergence acceleration algorithms [15] using 
singularity confinement. The integrability of the nonautonomous form of (2) (and its relation 
to cylindrical KdV) has been discussed by Nagai and Satsuma [16] in the framework of the 
bilinear formalism. 

We must point out here that the kind of initial conditions we choose, while influencing the 
speciflc degrees obtained, do not modify the conclusions on the type of growth. Let us 
illustrate this by choosing for (2) a staircase type of initial conditions where x~'^ — Pn/q, 
x]^'^ = Tnl q with the same convention as to the degrees of q and the p's, r's (but without 
the now unnecessary constraint pq = Tq). We flnd the degrees: 



11 16 
7 11 



1 1 



5 



where the underHned 1 corresponds to the origin. The growth is again quadratic and depends 
only on the sum N = n + m of the coordinates: = 1 + N{N — l)/2. 

Two more well-known discrete equations can be treated along the same lines. In the case of 
the lattice mKdV [14]: 

we obtain for constant z the same degree growth, = mn + 1, as for the potential lattice 
KdV. If we assume now a generic z we find the degrees: 



4 13 32 65 



13 21 31 



The degrees obey the recursion d"^^^ = c?^^^ + d^^+i + — ^ leading to an exponential 
growth with asymptotic ration (l+y/2). Requiring the degree of a;| to be 5 instead of 7 we 
find the condition 

m+l m _ m+1 m _ pi (r\ 

with solution z"^ = f{n)g{m). This condition is sufficient for the degrees of the nonau- 
tonomous case to coincide with those of the autonomous one. It is also precisely the one 
obtained in [15] using the singularity confinement condition. We believe that the Nagai- 
Satsuma approach [16] for the construction of double Casorati determinant solutions can be 
extended to the case of the nonautonomous lattice modified-KdV. 

The discrete sine-Gordon equation [17,18]: 
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(6) 



6 



in the autonomous case where z is a. constant leads to the degree pattern: 
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It can be represented by = mn + min(m, n) + 1. In the nonautonomous case of generic z 
we obtain the sequence of degrees: 
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obeying the relation dl^^^ = + + — (1 — 5^) leading again to exponential growth. 
The condition for a growth identical to that of the autonomous case is the same as (5). Thus 
equation (6) introduces a nonautonomous extension of the lattice sine-Gordon equation. 
(We intend to return to a study of its properties in some future work). We must point out 
that in the continuous limit, this nonautonomous form goes over to Wx,t = f {x)g(t) sin w. 
This explicit x and t dependence can be absorbed through a redefinition of the independent 
variables leading to the standard, autonomous, sine-Gordon, but no such gauge exists in the 
discrete case. 

We now turn to two discrete equations which are particular in the sense that they are not 
integrable through 1ST techniques but rather through direct linearisation. The first is the 
discrete Liouville equation [19]: 
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If we assume that z is a. constant we obtain the following degree pattern: 



By inspection we find d^^ = m + n. This result is not at all astonishing. As we have shown 
in [9], the degree growth of linearisable mappings is slower than that of the 1ST integrable 
ones. The same feature appears again here. The deautonomisation of (7) can proceed along 
the same lines as previously. For generic z™, the degrees are organised in a Pascal triangle 
and thus the growth is exponential. The condition for the growth to be identical to that of 
the autonomous case is again (5) and thus zl^ = f{n)g(rn). However, this nonautonomous 
extension is trivial: it can be obsorbed through a simple gauge transformation. Indeed, 
putting X = (j)X where = a{n)/3{in) with /(n) = a{n)a{n + 1), g{in) — P{m)P{m + 1) we 
can reduce equation (7) to one where z = 1. 

Finally, we analyse the discrete Burgers equation [19]: 

1 _|_ ytn m 

^-+i = a;-i±f!L^. (8) 

When z is a constant we find d'^ = m + (Notice that contrary to all the previous examples, 
in the case of Burgers equation m and n do not play the same role and thus a that is 
not symmetric in m, n is not surprising). For a generic z^, we find d^ = 2™", a manifestly 
exponential growth. The condition for the degree to grow like m + 1 is just 

C+i - C = (9) 

i.e. z^ = g{m). This leads to a nonautonomous extension of the lattice Burgers equation. 
Moreover this extension cannot be removed by a gauge. On the other hand, this extension 
is perfectly compatible with linearisability. Indeed, putting x'!^ = X^^^/X^ we can reduce 
it to the linear equation: 

X-+1 = /(m)(X- + ^(m)X- i) (10) 

where / is arbitrary and can be taken equal to unity. This nonautonomous extension is just 
a special case of the more general discrete Burgers: 

1 _i_ rurnrym ^ ' 

^ ^ In '^n 
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which can be Unearised through x"^ = ^^^X^^+^/X;^ to X;^+i = Vn (^^ + In^n^n+i) 
provided = 0:^7^1 and ct, and V' are related through q;^V^0^ = '0^+10^"'""'^ • We 
must point out here that the continuous Burgers equation also does possess a nonautonomous 
extension. It is straightforward to show that if is a solution of the equation 0^ = (f)'^(j)xx then 
the nonautonomous Burgers ut = + 2(f)uUx can be linearised to vt = (p^Vxx through 

the Cole-Hopf transformation u = (j)Vx/v. 

In this paper, we have applied the method of the slow degree growth to the study of the 
integrability of partial diflFerence equations. Our study has focused on well-known integrable 
lattice equations for which we have tried to provide nonautonomous forms. We have shown 
that using degree-growth methods it is possible to obtain integrable nonautonomous forms 
for most of the equations studied, and confirmed results previously obtained through the 
singularity confinement method. In the case of linearisable lattice equations, our results are 
the logical generalisation of the ones obtained for 1-dimensional mappings: the linearisable 
mappings have a degree growth that is slower than the one of the IST-integrable discrete 
equations. Our estimate of the degree growth was based on the direct computation of the 
degree for successive iterations and obtaining a fit of some analytical expression confirmed 
by subsequent iterations. It would be interesting, of course, to provide a rigorous proof 
of the degree growth following, for instance, the methods of [20]. However, this has not yet 
been carried through even for one- dimensional, nonautonomous mappings that are integrable 
through spectral methods. On the other hand, the proof of the degree growth for the cases 
where the equations are linearisable looks more tractable and we intend to address this 
question for both the one-and two-dimensional cases in some future work. 

The fact that we were able, through the adequate choice of initial data, to perform these 
calculations without being overwhelmed by their size is an indication of the usefulness of our 
approach. The study of degree growth, perhaps coupled with singularity confinement in the 
dual strategy we sketched in [10], can be a precious tool for the detection of integrability of 
multidimensional discrete systems. The interest of this method is not only that it can be 
used as a detector of new integrable lattice systems but also that it can furnish an indication 
as to the precise method of their integration. 
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